Stability of dark solitons in three dimensional dipolar Bose-Einstein condensates 
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The dynamical stability of dark solitons in dipolar Bose-Einstein condensates is studied. For stan- 
dard short-range interacting condensates dark solitons are unstable against transverse excitations 
in two and three dimensions. On the contrary, due to its non local character, the dipolar interaction 
allows for stable 3D stationary dark solitons, opening a qualitatively novel scenario in nonlinear 
atom optics. We discuss in detail the conditions to achieve this stability, which demand the use of 
an additional optical lattice, and the stability regimes. 



The physics of Bose-Einstein condensates (BECs) is, 
due to the interatomic interactions, inherently nonlinear, 
closely resembling the physics of other nonlinear systems, 
and in particular non linear optics. Non linear atom op- 
tics [H has indeed attracted a major attention in the last 
years, including phenomena like four- wave mixing [2| and 
condensate collapse [3|. One of the major consequences 
of nonlinearity is the possibility of achieving solitons, i.e. 
non-dispersive self-bound solutions, in quasi-one dimen- 
sional BECs. Bright solitons have been reported in BECs 
with negative s-wave scattering length a < (the equiva- 
lent of self- focusing nonlinearity) [4j . Dark solitons (DSs) 
have been realized as well in BECs with a > (self- 
defocusing nonlinearity) In addition, optical lattices 
have allowed for the observation of gap solitons Q . 

The stability of solitons depends crucially on the quasi- 
one dimensionality of the systems, which for the case of 
BEC solitons demands a sufficiently strong transversal 
confinement of the condensates Q. For the particular 
case of dark solitons, if the transversal size of the sys- 
tem becomes comparable to the width of the dark soliton 
(typically provided by the healing length of the system, as 
discussed below), then the dark soliton plane becomes dy- 
namically unstable. This dynamical instability (so-called 
snake instability) has been previously studied in the con- 
text of non linear optics In the context of BEC, it has 
been shown that this instability leads to a strong bending 
of the nodal plane, which breaks down into vortex rings 
and sound excitations 0], as experimentally observed in 
Ref. 0. 

Nonlinear phenomena constitute an excellent exam- 
ple of the crucial role played by interactions in quantum 
gases. Until recently, typical experiments involved parti- 
cles interacting dominantly via short-range isotropic po- 
tentials, which, due to the very low energies involved, are 
fully determined by the corresponding s-wave scattering 
length. However, recent experiments on cold molecules 
Rydberg atoms 12 1, and atoms with large mag- 



netic moment [13( , open a fascinating new research area, 
namely that of dipolar gases, for which the dipole-dipole 
interaction (DDI) plays a significant or even dominant 



role. The DDI is long-range and anisotropic (partially at- 
tractive) , and leads to fundamentally new physics in ultra 
cold gases (T3 | . Time-of- flight experiments in Chromium 
have allowed for the first observation of DDI effects in 
cold gases jlH, which have been remarkably enhanced 
recently by means of Feshbach resonances 

Dipolar gases present a rich non-linear physics, since 
the DDI leads to non-local non-linearity, similar as that 
encountered in plasmas fl7f. nematic liquid crystals [l8| , 
thermo-optical materials [19j | and photorrefractive crys- 
tals [13] ■ Nonlocality leads to a wealth of novel phenom- 
ena in nonlinear physics, as the modification of modula- 
tion instability |21| , the change of the soliton interaction 
[22T |. and the stabilization of azimuthons [23}. Particu- 
larly interesting is the possibility of stabilization of local- 
ized waves in cubic nonlinear materials with a symmetric 
nonlocal nonlinear response [HI]- Multidimensional soli- 
tons have been experimentally observed in nematic liq- 
uid crystals 25[ and in photorrefractive screening solitons 
[2^ . Recently we showed that under realistic conditions, 
2D bright solitons may be generated in dipolar BEC 27 1 . 
Interestingly, due to the non-local nature of the DDI, the 
soliton-soliton scattering presents strong inelastic reso- 
nances, as well as spiraling motion [2SL similar to that 
observed in photorrefractive materials 26j . 

In this letter, we show that the long-range charac- 
ter of the DDI may have striking consequences for the 
stability of dark solitons in dipolar BECs. Contrary to 
usual BECs, for which, as mentioned above, dark nodal 
planes become unstable when departing from the one- 
dimensional condition, the DDI may stabilize dark nodal 
planes even if the transversal size of the condensate be- 
comes arbitrarily much wider than the condensate heal- 
ing length. This stabilization is purely due to the long- 
range character of the DDI. We study in detail the condi- 
tions for this stabilization, and the stabilization regimes. 

In the following, we consider a dipolar BEC of parti- 
cles with mass m and electric dipole d (the results are 
equally valid for magnetic dipoles) oriented in the z- 
direction by a sufficiently large external field, and that 
hence interact via a dipole-dipole potential: Vd(r) = 



2 



ad 2 (I — 3 cos 2 (#))/r 3 , where 9 is the angle formed by 
the vector joining the interacting particles and the dipole 
interaction. The coefficient a can be tuned within the 
range — l/2<a<lby rotating the external field that 
orients the dipoles much faster than any other relevant 
time scale in the system 29] . At sufficiently low temper- 
atures the physics of the dipolar BEC is provided by a 
non-local non-linear Schrodinger equation (NLSE) of the 
form: 



d 

dt y ' 



2m 



V ol {x,y)+g\*{?,t)f 



dr'V d (f-f')\^^,t)\ 2 



(1) 



where g — 4irh 2 a/m, with a the s-wave scattering length 
(we consider a > 0) and m the particle mass. For rea- 
sons that will become clear below, the BEC is assumed 
to be in a 2D optical lattice, V i(x,y) — sEn(sin 2 (qix) + 
sin 2 (qiy)), where Er = h 2 qf/2m is the recoil energy, qi is 
the laser wave vector and s is a dimcnsionless parameter 
providing the lattice depth. In the tight-binding regime 
(i.e. for a sufficiently strong lattice but still maintaining 
coherence), we may write ^(r, t) = J^ijfij (x, y)ipi.j {z, t), 
where fij(x, y) is the Wannier function associated to the 
lowest energy band and the site located at (bi,bj), with 
b = ir/qi- Substituting this ansatz in Eq. ([1]) we obtain 
a discrete NLSE [3(J. We may then return to a contin- 
uous equation, where the lattice is taken into account in 
an effective mass along the lattice directions and in the 
renormalization of the coupling constant 31]. 



d 

dt K ' ' 



h 2 „ 2 h 2 d 2 _ 1T , W . l2 



dfV d (f-r')\^(r',t)f 



(2) 



where g = b 2 g J f(x,y) 4 dxdy + g d C j3g], with g d = 



a8nd 2 /3, 



h /2b J is the effective mass, and 



J = J dxdyf i j(x,y)[~(h 2 /2m)S7 2 L + V ol {x,y)}f Vj >{x,y), 
for neighboring sites and The validity of 

Eq. ([2]) is limited to radial momenta k p <C 2w/b, in which 
one can ignore the discreteness of lattice. In the following 
we use the convenient dimensionless parameter (3 — g d /g, 
that characterizes the strength of the DDI compared to 
the short range interaction. The Fourier transform of the 
DDI: V d {k) = .gd[3cos 2 6 k - l]/2, with cos 2 9 k = fc 2 /|fc| 2 , 
is needed later for the calculations. 

Due to its partially attractive character, the stability of 
a dipolar BEC is a matter of obvious concern [3] . A Bo- 
goliubov analysis of an homogeneous dipolar condensate 
reveals that the dispersion relation for quasiparticles is 
of thejbrm e(k) = [E kin (k)[E kin (k) + E^k)]] 1 / 2 , where 
Ehinik) = h 2 k 2 /2m* + H 2 k 2 /2m is the kinetic energy, 
and Ei nt (k) = 2(g + V d (k))no is the interaction energy, 
which includes both short-range and dipolar parts. Note 



that V d (k) may be positive or negative, and hence for low 
momenta (phonon excitations) the dynamical instability 
(phonon instability) is just prevented if — 1 < (3 < 2. If 
g d > 0, phonons with k lying on the xy plane are unsta- 
ble if f3 > 2, while f or g d < phonons with k along z are 
unstable if j3 < — 1 [33j. 

In this paper, we are particularly concerned about the 
stability of a dark soliton in a three dimensional dipo- 
lar BEC. We assume that the dark soliton lies on the xy 
plane, hence the solution can be written in the following 
form: ^o(f, t) = ipo(z) expj— ifit/h], where [i is the chem- 
ical potential. Introducing this expression into Eq. ([1]) 
we obtain a one-dimensional NLSE in z of the form: 



'^ 2+9lMz)l 



Mz)- (3) 



Note that, since ipo is independent of x and y, in this 
equation the DDI interaction just regularizes the value 
of the local coupling constant g = g + g d . This equation 
allows for a simple solution describing a dark soliton, 
■0o = v / ^o^ ann ( z /0i where ( = h/^/mgno is the corre- 
sponding healing length and uq is the bulk density. Note 
that, interestingly, due to the modification of the local 
coupling constant the size £ of the dark soliton depends 
on the DDI. 

We are especially interested in the dynamical stability 
of these nodal planes. To study this stability we perform 
a Bogoliubov analysis, considering a transversal pertur- 
bation of the dark soliton planes: \&(r, t) = ^o(r,t) + 
x(r, t) exp(— ifjit/K), where x(?,t) = u{z) exp[z(gx — 
ct/ti)] + v(z)* exp[— i(qx — e*t/H)], where q is the momen- 
tum of the transverse modes with energy e. Introducing 
this ansatz into Eq.|T]) and linearizing in x, one obtains 
the corresponding Bogoliubov-de Gennes (BdG) equa- 
tions for the excitation energies e and the corresponding 
eigenfunctions f± = u ± v. 



£/-(*) = 



2m \ dz 2 



q - fi + 3gip (zY 



/+(*) 



--g d qMz) I dz' exp(-g|z - z'\)ip (z')f+(z'), (4) 



tf_ f d 2 

2m \ dz 2 



Note that the dipole interaction has two main effects. 
On one side, as mentioned above, it leads to a regularized 
g. On the other side, it introduces a qualitatively new 
term in the second line of Eq. 2J Whereas the first ef- 
fect just leads to a quantitative modification of the dark 
soliton width, the second effect is a purely dipolc-induced 
non local effect, which, as we show below, may lead to re- 
markable consequences for the dark soliton stability. For 
every transverse momentum q we determine the lowest 
eigenenergy, which provides the dispersion law e(q). 
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Figure 1: Imaginary part of the excitation energies of a DS 
for m/m* — 1, and ,0 = (triangles), —0.5 (squares) and 
1 (circles). Solid lines correspond to the analytical result for 
low momenta. The real part of the excitation energies is equal 
to zero for the range of momenta considered in the figure. 

When (3 = (no DDI) and m/m* = 1 (no lattice), 
we recover the same BdG equations discussed in the con- 
text of standard short-range interacting condensates Q • 
It has been shown that in that case the dispersion law 
e(q) is purely imaginary for qC, < 1 (Fig.[T]). Hence dark- 
soliton planes in homogeneous three-dimensional short- 
range interacting condensates, are dynamically unstable 
against transverse modulations. This so-called snake in- 
stability has been experimentally observed in non linear 
optics |8[ and recently in the context of BEC [lfj . In the 
latter case, the dark-soliton plane bends and decays into 
more stable structures such as vortex rings and sound ex- 
citations [9( . The stabilization of the soliton demands a 
strong radial confinement i-e. the dark-solitons can- 
not be considered any longer as three-dimensional, but 
on the contrary acquire a one-dimensional character. 

In the presence of DDI ((3 ^ 0) but without lattice 
(m/m* = 1), the transverse instability persists, since e(q) 
remains purely imaginary for q( < 1. For (3 > ((3 < 0), 
\e(q)\ decreases (increases) when \(3\ grows (Fig. [I]). The 
situation changes completely when m/m* is sufficiently 
small since, as discussed below, e(q) is purely real and, 
as a consequence, the dark soliton becomes dynamically 
stable (see FigJ3]). As we detail in the following, this re- 
markable fact can be clearly understood by means of a 
transparent physical picture. First we notice that for low 
momenta the spectrum is always linear in q, suggesting 
the idea that for low momenta the system may be de- 
scribed by an elastic model of the two dimensional nodal 
plane of the dark soliton. The Lagrangian density for the 
nodal plane reads 

OK(§)'-fw (o» 

where <f> is the displacement field of the nodal plane from 
the ground state, M is the mass for unit area and a 
plays the role of a surface tension. The mass M of the 
soliton can be easily calculated expanding the energy 



Figure 2: Analytical (solid line) and numerical (empty circles) 
results for the stable/unstable phonon regimes as a function 
of two paremeters (3 and m/m* . 



of a moving soliton up to second order in the velocity. 
We obtain M = —4hno/c, where c = \J gn^/m is the 
sound velocity. Notice that M < since the dark soli- 
ton represents an absence of atoms. The surface tension 
can be calculated inserting a suitable variational ansatz 
^via(r) = Y / n u 'tanh((z — \/2a cos(qx)) / X) (describing a 
transverse modulation of the nodal plane with amplitude 
a and momentum q) in the energy functional and ex- 
panding up to second order in a and q: 



4n h 2 
3(m* 



(7) 



This expression can be considered as one of the main 
results of this Letter. The nature of the eigenmodcs 
uj 2 = (a/M)q 2 crucially depends on the sign of cr/M. 
In the absence of DDI ((3 = 0) , a is always positive and 
hence the modes are purely imaginary. The dark soli- 
ton is dynamically unstable against the above mentioned 
snake instability. Note that for (3 = and m/m* = 1 
our result coincides with the one found in Ref . [7|] . In the 
absence of an additional optical lattice, the dynamical 
instability of the DS at low q dissapears for (3 > 2, i.e. 
for situations for which the homogeneous dipolar BEC as 
a whole is itself, as commented above, unstable against 
local collapses. Increasing the depth of the lattice poten- 
tial reduces the role of the kinetic energy term (m/m*)q 2 
in Eqs. ^ and ([5]) (or equivalently reduces the first term 
in Eq. ([7])) and hence enhances the role of the DDI. A 
sufficiently large DDI or small m/m* such that 



m 

TO* 



< 



3/'? 



2(1 + 0) 



(8) 



leads to stable low-energy phonons with q — > 0. We have 
evaluated from a direct numerical calculation of the BdG 
equations (|4]) and ([5]) the stability threshold at which 
e(q — > 0) becomes real. Fig. [2] compares our numerical 
and analytical results, which are in excellent agreement. 

When m/m* decreases further or (3 grows, a wider 
regime of low momenta is stabilized (Fig. [3]). Note in 
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Figure 3: Real part of the excitation energies of a DS for 
m/m* — .1, and (3 — 1.6. Solid line corresponds to the ana- 
lytical result for low momenta while empty circles correspond 
to numerical results. The imaginary part of the excitation 
energies is equal to zero for the range of momenta considered 
in the figure. 



Fig. [3] that the dispersion law at low momenta is very 
accurately described by our analytical results. A suffi- 
ciently strong optical lattice and large DDI can stabi- 
lize all the modes with momenta up to q^/m/m*C ~ 
1. Although we observe instabilities for momenta 
q^Jm/m*C, ~ 1, this large- momentum instability is typi- 
cally irrelevant, since for sufficiently small m/m* it con- 
cerns momenta much larger than the lattice momentum. 
Although our effective mass theory breaks down for such 
momenta, it becomes clear that such high momentum 
instabilities are physically prevented by the zero point 
oscillations at each lattice site. 

Summarizing, contrary to short-range interacting 
BECs, where stable dark solitons demand a sufficiently 
strong transverse confinement, dipolar BECs allow for 
stable dark solitons of arbitrarily large transversal sizes. 
We have obtained the stability conditions, which demand 
a sufficiently large dipole and a sufficiently deep optical 
lattice in the nodal plane. We presented analytical re- 
sults for the lowest part of the spectrum which agree with 
those obtained numerically. The stabilization of nodal 
planes is purely linked to the long-range nature of the 
DDI, opening a qualitatively new scenario in non linear 
atom optics. 
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